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We discuss ILC measurements for a specific MSSM scenario with CP phases, 
where the lightest neutralino is a good candidate for dark matter, annihilating 
efficiently through t-channel exchange of light staus. These prospective (CP-even) 
ILC measurements are then used to fit the underlying model parameters. A collider 
prediction of the relic density of the neutralino from this fit gives 0.116 < VLh? < 0.19 
at 95% CL. CP-odd observables, while being a direct signal of CP violation, do not 
help in further constraining Qh?. The interplay with (in)direct detection of dark 
matter and with measurements of electric dipole moments is also discussed. Finally 
we comment on collider measurements at higher energies for refining the prediction 
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1 Introduction 

One of the prime motivations for a high-energy and high-luminosity e + e~ linear collider 
is the possibility to do precision measurements of new particles beyond the Standard 
Model. The lightest of these new particles is often stable by virtue of a new discrete 
symmetry and hence a candidate for the dark matter (DM) of the Universe. One therefore 
hopes to be able to precisely determine the properties of the DM in the laboratory, 
and in particular to make a "collider prediction" of its relic abundance, which can be 
tested against cosmological models. For such a collider prediction to be of interest, it 
must be as precise as the value obtained from cosmological observations. This means 
a precision of about 10% (for 95% CL) to match WMAP+SDSS [1-3] or few percent 
to match expectations at the PLANCK satellite [4]. Moreover, if the annihilation cross 
section of the DM candidate is known, one can also predict direct and indirect DM 
detection cross sections as functions of astrophysical quantities such as galactic densities, 
see e.g. [5, 6]. 

The possibility to make collider predictions of the cross sections for annihilation of dark 
matter candidates has been examined within specific supersymmetric scenarios. Within 
the constrained minimal supersymmetric standard model (CMSSM), which has only a 
handful of free parameters, it has been shown in a case study [7] that the LHC could 
make a prediction for the DM relic density of the order of 10%, assuming the standard 
cosmological scenario; the ILC could reach a precision of few percent [8] in several scenar- 
ios. The particular case of stau co-annihilation in the CMSSM was investigated for the 
LHC in [9] and for the ILC in [10] with similar conclusions. In the general MSSM, it was 
shown that the LHC might match roughly the WMAP+SDSS precision in a favourable 
scenario [11] while the ILC could achieve much better precision [12]. These conclusions, 
however, depend very strongly on the scenario considered; many remain challenging for 
both the LHC and the ILC [12], see also [13]. Moreover, the studies mentioned here mainly 
concentrated on a few scenarios that provide the correct amount of neutralino annihila- 
tion, consistent with the WMAP+SDSS range. Furthermore -and more importantly- 
they assumed that CP is conserved, although CP-violating (CPV) phases are generic in 
the MSSM. 

It is well known that CPV MSSM phases can have an important effect on the neutralino 
annihilation cross sections [14-18]. The consequences of phases for direct and indirect 
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detection were examined in [15,16,19-21]. In [18], we performed a comprehensive analysis 
of the impact of CP phases on the relic density of neutralino dark matter, taking into 
account consistently phases in all (co-) annihilation channels, and carefully disentangling 
CPV effects due to modifications in the couplings from pure kinematic effects. We found 
variations in Qh 2 solely from modifications in the couplings of up to an order of magnitude. 
We concluded that the determination of the relevant couplings (including CPV phases!) 
can be as important for the prediction of Qh 2 as precision measurements of masses, i.e. 
pure sparticle spectroscopy. 

In this paper, we therefore consider a particular scenario of the CPV MSSM, taken 
from [18], and investigate i) which measurements are possible at the ILC, ii) to which 
precision the underlying MSSM parameters and the neutralino relic density can be inferred 
from these measurements, and iii) what are the implications for direct and indirect DM 
detection and measurements of electric dipole moments (EDMs). 

Our scenario has light gauginos and staus with masses below 200 GeV. The lightest 
super symmetric particle (LSP) is the lightest neutralino with a mass of ~80 GeV; it is 
dominantly bino. The two staus have masses of about 100 GeV and 180 GeV, with a 
strong mixing between the left- and right-chiral states. The neutralino LSP annihilates 
predominantly into tau pairs, with the annihilation cross section being sensitive to the 
stau mixing. Although the f\ is very light, the scenario does not rely on coannihilation 
but on t-channel stau exchange. We therefore refer to it as the "stau-bulk" scenario. Such 
a scenario occurs in the CP-conserving MSSM only for M\ < 0. 

We further impose that the sfermions of the first and second generation are heavy in 
order to avoid the strong EDM constraints. The resulting mass pattern, light staus but 
TeV-scale selectrons and smuons, is not found in SUSY models where universality among 
scalar masses is imposed at a high scale. It is important to note that our scenario, despite 
having several particles below 200 GeV, is quite challenging for colliders. At the LHC, 
SUSY events are dominated by squark and gluino production followed by cascade decays 
leading to jets plus r's plus E™ 1SS - At the ILC, two neutralinos, the lighter chargino, and 
the two staus are within kinematic reach with y^i = 500 GeV. However, production of 
these sparticles again only leads to taus plus missing energy. 

The paper is organised as follows. Section 2 describes the detailed setup of our "stau- 
bulk" scenario. Expectations for ILC measurements are given in Section 3. In Section 4, 
we present our results concerning the determination of the model parameters and the 
prediction for the relic density of dark matter. In Section 5 we then discuss CP-sensitive 
observables, and in Section 6 additional possibilities to test and constrain the model at 
higher energies, through EDM measurements or through direct DM detection experiments. 
Finally, Section 7 contains our conclusions. 

2 Setup of the stau-bulk scenario 

In the MSSM, the parameters that can have CP phases are the gaugino and Higgsino mass 
parameters, Mj = \M i \e l ^ i , with i = 1,2,3, and fi = |/i|e"^ M , and the trilinear sfermion- 
Higgs couplings, Af — \Af\e t<l,f . Not all of these phases are, however, physical. The 
physical combinations are Arg(Mj/x) and Arg(A//i). Allowing for CP- violating phases, in 
particular non-vanishing 0/, induces a mixing between the two CP-even states h°, H° and 
the CP-odd state A . The resulting mass eigenstates hi, h 2 , h 3 (m hl < m h2 < m hs ) are 
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Table 1: Masses of particles for the input parameters of eq. ([T]). 



no longer eigenstates of CP. Therefore the charged Higgs mass, , is typically used as 
input parameter in the CPV-MSSM. 

In this paper we investigate the "stau bulk" region of [18], which appears for light 
staus and large phase of Mi. The point is the following. In the conventional case, there 
are mainly two mechanisms that can make a light bino-LSP annihilate efficiently enough: 
resonant annihilation through the light Higgs, or co-annihilation with a sparticle that 
is close in mass -in the CMSSM typically the lighter stau. The so-called "bulk region" 
where the LSP annihilates through t-channel exchange of very light sleptons is largely 
excluded by LEP. For large phases of M\ we have found, however, that the couplings of 
the neutralino to staus can be sufficiently enhanced such that a new region opens up, 
where the Xi~^i mass difference is too large for co-annihilation to be efficient but the LSP 
annihilates into taus through t-channel exchange of both f\ and f%. 

To define a benchmark point for such a scenario, we choose the following input pa- 
rameters at the electroweak scale: 

Mi = 80.47 GeV, M 2 = 170.35 GeV, M 3 = 700 GeV, fa = 180°, 

M u = 138.7 GeV, = 135.2 GeV, A T = 60 GeV, <p T = 0, (1) 

fi = 600 GeV, tan (3 = 10, </v = 0. 

All other parameters, i.e. sfermion masses, A t> b, and m#+ are set to 1 TeV for simplicity. 
This way EDM constraints are avoided when varying </>i and A T ; all other phases are 
assumed to vanish. 

We use CPsuperH [22] to compute Higgs and sparticle masses and mixing angles, and 
micrOMEGAs 2 . 1 [23,24] to compute the relic density, EDMs, and (in)direct detection cross 
sections. The mixing angle in the stau sector writes 

tan 2 0, = -^r(A T -^n(3) 
mi — mi 

and for our scenario is completely dominated by the /it an/? term. The mass spectrum 
resulting from eq. (CQ) is given in Table [U The particles accessible at ILC with = 
500 GeV are Xii X°> Xii ^2 and v T . Production cross sections and branching ratios are 
discussed in Section 3. At this stage we just note that the scenario is rather challenging 
to resolve experimentally because SUSY events involve only r's and missing energy. 

The relic density of the Xi is ^h 2 = 0.130 at the nominal point eq. ([T]). As mentioned 
above, the dominant annihilation channel is into tau pairs (more than 95% of the total 
contribution) through t-channel exchange of staus. The contribution of both fi and f 2 
is crucial in bringing Qh 2 to the desired range, 0.094 < Qh 2 < 0.136 [25]. Indeed, for 
0i = (or Mi positive) one would have Qh 2 = 0.167. Note also that the mass splitting 
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Figure 1: Dependence of Qh 2 on (dashed red line) and m^j (full blue line) for the 
nominal point eq. (OQ). The horizontal lines indicate 5Q/Q — ±5%. 

between the stau-NLSP and the LSP is too large for any significant contribution from 
co-annihilation. 

The precision with which Qh 2 can be inferred from ILC measurements therefore de- 
pends not only on the accuracy of the sparticle spectroscopy (mass measurements) but on 
the determination of all parameters of the neutralino sector (Mi, M2, //, tan/5, <pi), which 
determine the LSP mass and couplings, and the four parameters of the stau sector (Mf L , 
Mf R , A T , (f) T ). The dependence on <pi, and to a much lesser extent <fr T , originates from the 
Xifi 2 T couplings. In addition, particles which are too heavy to be produced at ILC could 
have some influence. Indeed, the exclusion of selectrons and smuons up the the kinematic 
limit, mj > 250 GeV, leaves an uncertainty in Qh 2 of about 7%. The influence of heavy 
Higgs states gives 5Q/Q ~ 5% for mh 23 — 250 GeV. The combined effect from sleptons 
and Higgses is actually smaller, because their contributions work against each other, so 
we have <5f2/f2 ~ 7% as overall uncertainty from the unknown part of the spectrum^] The 
dependence on and m h2Z is illustrated in Fig. [TJ 

We assume here that the mass scale of the squarks and gluino is known from LHC. 
The LHC can also provide additional information on the heavy Higgs states: a charged 
Higgs with mass up to ~ 500 GeV could be discovered/excluded at the LHC in the 
channel gb — > H~t with H~ — > tv for tan/3 = 10 [26]. If the neutral heavy Higgses weigh 
more than 250 GeV, the LHC has limited discovery potential using SM decay modes at 
intermediate values of tan/3. However, a recent analysis [27] has shown how to exploit 
the Higgs into chargino/neutralino decay modes into four leptons (electrons or muons) in 
this region. Although no dedicated analysis has been performed for decays into taus, a 
similar coverage as in [27] can be expected in our scenario [28]. 

1 Notice, however, that the measurement of the total e + e~ — ► t + t~ cross section above the Xi~Xr 
threshold, being dominated by chargino-pair production, provides an indirect constraint on the v e mass, 
and hence also on that of e^, of m^ c ^ L ^ 900 GeV. 
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3 ILC measurements 



3.1 Event generation and analysis 

The analysis presented in this section assumes a detector as described in the Tesla 
TDR [29]. The simulations are based on Simdet 4.02 [30] with an acceptance 9 > 
125 mrad and e, 7 veto of 9 > 4.6 mrad. Background and signal events are generated with 
Pythia 6.2 [31] including beam polarisations {V e -,V e +) = (0.8,0.6), QED radiation and 
beamstrahlung a la ClRCE [32]. r decays and polarisation are treated with Tauola [33]. 
The SM background includes e + e _ — > W + W~ and the negligible pair production e + e~ — > 
r + T~(-f). The 77 background is also negligible with a(e + e~ — > e + e~r + r~) = 4.5 ■ 10 5 fb 
and an acceptance < fewlO" 6 . 

To measure the particle masses, two methods are available: threshold scans and end- 
point methods. Both will be useful in our case. In slepton decays to LSP, t~ — > t~ Xi, 
the endpoints, i.e. maxima and minima E+/- in the flat lepton energy spectra: 



where (3 = * 1 — 4m~/s, can be used to extract the slepton and neutralino mass. 



The mixing angle in the stau sector, eq. ([2D, can either be determined from measure- 
ments of the polarized cross section or from the measurement of the r polarisation V T in 
the decay f\ — > TXi [34-36] . Note that V T depends not only on the mixing anlge 9f of the 
staus but also on the mixing of the lightest neutralino. 

The selection criteria used in e + e~ — > f^fi — > T + T~E mms are essentially to require two 
acoplanar r jets (A0 < 180° and m T < 2 GeV) being produced centrally (|cos#| < 0.75). 
In addition, topology dependent p± cuts are applied to reject 77 reactions [10]. The 
overall efficiency for fi-pair reconstruction is e ~ 0.32 at ^/s = 280 GeV and varies only 
slowly with the cms energy. The same selection criteria are applied to detect the other 
sparticle decays during an energy scan. 

Despite the large statistics, the leptonic r decays are not useful due to the large 
WW background. We therefore only analyse the decay modes r — > 7iv T , pu T1 ?>ixv T with 
branching ratios B(r -> ttz/ t ) = 0.111, B(r pu T ) = 0.254, B(r -> 3ttz/ t ) = 0.194. 

3.2 Masses from threshold scan 

At = 500 GeV, all SUSY signals are completely dominated by the TrE miss topology. 
The challenge is then to disentangle the various sources. The strategy we propose is to 
scan downwards in energy in steps of 10 GeV in order to find thresholds. This is done 
using different beam polarisations e^ej [crl] and e^e^ \o~lr] where V e - = ±0.8 and 
V e + = =f0.6. The f^f~ pair production is best detected in the orx polarisation mode, 
while cflr is dominated by XiXi production The results of a scan of the visible cross 
sections e + e~ — > rrE miss are shown in Fig. [21 Clearly, integrated luminosities of 2 fb -1 
per step in energy are sufficient. Fitting the excitation curves, Off ~ j3 3 of Fig. [2^ and 
a x + x~ ~ °f Fig- Eb, allows a determination of the sparticle masses. Note that the 
mixed f\fi production is hardly detectable and can only be accomodated in a global fit 
including all channels. The results of a fit for f 1? xt an d T2 masses together with the total 
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Figure 2: Visible cross section as a function of cms energy of the inclusive reaction e^e^ — > 
t + t~ E miss (left) and e^e^ — > t + t~ E miss (right). The red (blue) points correspond to 
integrated luminosities of C = 2fb~ (10fb~ ), the contributions of W + W~ (dash-dot), 
fi?i (blue), f x T2 (pink), XiXi (green) and t^tT" are indicated. 

integrated luminosity used in each case are listed in Table El The expected accuracies can 
easily be scaled with the luminosity. 

reaction m [GeV] C [fb _1 ] 

e l e R ~ Wi 100.92 ±0.40 50 

e R e l^XxXx 164.88 ±0.015 30 
e\e- R ^T2f 2 188.2 ±9.1 100 

Table 2: Results of mass fits from the excitation curves indicating the reactions, beam 
polarisation and assumed integrated luminosity 



3.3 LSP mass from stau decays 

An alternative method for measuring the f\ mass and/or the Xi LSP mass is to use the 
energy spectra in f\ — > t\\ decays. However, the r energy spectrum reconstructed from 
the r decay products is no longer flat due to undetected neutrinos, and the extraction of 
the endpoint energies, see eq. ©, gets more involved [35,36]. The upper endpoint energy 
E + can be identified with the maximum energy of the decay products and can still be 
determined rather well. The lower endpoint energy EL is in general completely distorted. 
It may be reconstructed from the pion energy E n spectrum of r — > ix u T , but the expected 
precision is poor due the low branching ratio and a polarisation V T dependent shape of 
the energy distribution. Therefore, determining both stau and neutralino masses from 
the energy spectra alone is not very meaningful and we use the additional information on 
from the threshold scan. 

In order to study the properties of stau decay spectra, the maximal energy for f^ff 
production at y/s = 280 GeV just below the threshold for other sparticles is chosen. A fit 
to the E w spectrum (not shown) yields E_ = 7.8 ± 0.50 GeV and E + = 42.8 ± 0.30 GeV. 
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Figure 3: Energy spectrum Ej et from r 

including background from W + W~ — > t + t & m i SS ] ^ne aata points represent a simulation 
assuming C = 200 fb" 1 at v/i = 280 GeV, the blue histo gram corresponds to a very high 
statistics sample, the green curve is a fit to an analytic formula (left). Correlation m^o 
vs rrtfj using the stau mass from the excitation curve and the end point energy from the 
Ej et spectrum (right). 



p u T + 3tt u t of the reaction e J [e R — > f^f x 



In r — > p v T — > 7T ± 7r° v T and r — > 37r v t — > 7r ± 7r ± 7r =F z/ r + 7r ± 7r 7r° v T the energies E p and 
-E^ have the advantage of being independent of V T . The analysis of Ej et (= E p + £"3^) 
spectrum is shown in Fig. [31 a fit to the energy distribution using an analytic formula 
gives the upper endpoint E + = 42.8 ± 0.20 GeV. Applying eq. together with the 
results from the excitation curve leads to a LSP mass of m^o = 80.67 ± 0.35 GeV and a 
stau mass of = 100.92 ± 0.39 GeV. The m^o — correlation is shown in Fig. [31 
Obviously an improved accuracy on will result in a smaller error on the LSP mass. 



3.4 Stau mixing angle 

The stau mixing angle can be determined from measurements of polarized cross sections. 
Again, the energy for f\ pair production is chosen below the threshold for other sparticles 
production, that is y/s = 280 GeV. At this energy we expect ctrl = 155.5 ± 2.2 fb where 
the error corresponds to an integrated luminosity of C = 200 fb -1 and includes an esti- 
mate of systematics from acceptance calculations, branching ratios and WW background 
subtraction; for comparison see the visible cross section in Fig. [2j 

The polarised cross section <jrl dependence on the mixing angle cos 29 f [36] is shown in 
Fig. HI Fitting both the mixing angle cos 29 f and the mass to the cross sections in the 
continuum and at threshold yields cos20 f = -0.065 ±0.028 and m fl = 100.92 ±0.40 GeV. 
The 771.^—008 2^ correlation contour is shown in Fig. HI The accuracy can be improved 
by adding another cross section measurement with different polarisation, e.g. <jlr = 
157.3 ± 3.5 fb assuming C = 200 fb -1 . Using all observables, (Jrl-, o~lr and mf^athr), one 
expects a precision of Sm fl = 0.35 GeV and 5cos29f = 0.017. 
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cos 20 T cos 20 T 

Figure 4: Dependence of fifi polarised cross section ctrl on the f mixing angle cos 26 f 
at y^i = 280 GeV, the blue band indicates a measurement assuming C = 2001b" 1 (left): 
correlation m fl vs cos 26 f using the stau mass from the excitaion curve the cross section 
a RL (right). 

3.5 Tau polarisation 

The polarisation V T of the tau stemming from the fj — * rx? decay gives additional infor- 
mation on the stau and neutralino mixings [34]. Since cos 26 f is obtained from polarized 
cross sections as explained above, V T is in particular useful to constrain the gaugino- 
higgsino composition of the LSP. Here, we exploit the decay r — > pv T — > tt7t u t . The 
ratio z n = E n /E p is indeed sensitive to V T and independent of m^. For a right-handed 
r (V T = +1), the p is longitudinally polarized and the distribution da/dz n oc (2z n — l) 2 
is peaked near z n — > or 1, while for a left-handed r (V T = — 1), the p is transversely 
polarized and the distribution da/dz n oc 2z 7T (l — z n ) is peaked at z n = 0.5. The analysis 
of the z v spectrum is shown in Fig. El The results of a fit to the fi polarisation leads to 
V T = 0.64 ± 0.035 for an input value V T = 0.641. 

In summary, simulations of stau production e + e~ — > fjfj and chargino production 
e + e~ —> XiXi under realistic experimental conditions assuming (V e -,V e +) = (0.8,0.6) 
show that the stau, neutralino and chargino masses as well as the polarisation V fl _, T ^o and 
mixing parameter 6f can be accurately determined with moderate integrated luminosity. 
The possible ILC measurements are summarized in Table [3] 



channel 


observables 


f+ff 


m fl = 100.92 ± 0.40 GeV 


Al 


= 80.67 ±0.35 GeV 




cos 20 f = -0.065 ±0.028 


v T = 


0.64 ±0.035 




m f2 = 176.9 ±9.1 GeV 






xfxl 


m-± = 164.88 ± 0.015 GeV 







Table 3: Summary of achievable precisions of ILC measurements. 
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Figure 5: Energy spectrum z n = E n /E p from r — > p v T — > 7T7r° z/ T of the reaction 
e L e ij ~~ * ^i~^r including background from W + W~ — > r + z/ t _ za The data points represent 
a simulation assuming C = 200 fb" 1 at yfs = 280 GeV, the blue histo gram corresponds 
to a very high statistics sample, the blue curve is a fit to the data with a fi polarisation 
V T = 0.64 ±0.035. 

4 DM properties: fit to ILC observables 

Having established how precisely one could measure various CP-even observables at ILC, 
we now estimate how this would constrain the parameters of the underlying model, and 
the properties of the dark matter candidate. To this aim we perform a fit to the six 
observables listed in Table [3] assuming Gaussian errors. We also compute the total rr 
SUSY cross section at = 400 GeV with polarized beams and add it to the fit: ct(tt) = 
3.220±0.046 pb (for 10 fb" 1 ). Here we assume a systematic error as large as the statistical 
error. As free parameters we take 

M u n, tan/?, M Ls , M^, A T , 1; <f> T . (4) 

Owing to the extremely small experimental error on m-±, we do not include M 2 as an 
independent variable in the fit but compute it from /x and tan (3 such that it matches the 
measured value of m^± . The other SUSY parameters are fixed to the values specified in 
Sect. |2j We do not include the light Higgs mass in the fit because of the large parametric 
uncertainty from not knowing the parameters of the stop sector. This will be discussed 
in more detail after having summarized our main results. 

To probe efficiently the multidimensional parameter space, we perform a Markov 
Chain [37] Monte Carlo (MCMC) analysis using a Metropolis-Hastings algorithm [38-40]. 
This algorithm generates a candidate state x c from the present state x l using a proposal 
density Q(x';x c ). The candidate state is accepted to be the next state x t+1 in the chain 
if 

= P(a; c )g(x t ;x c ) 

where P(x) is the probability calculated for the state x, is greater then a uniform random 
number a — U(0, 1). If the candidate is not accepted, the present state x t is retained and a 
new candidate state is generated. For the proposal density, we use a Gaussian distribution 
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Figure 6: Results of the MCMC fit; in (a) x 2 as a function of /i, in (b) the la (black) and 
2a (green) allowed regions projected onto the \x versus tan (3 plane. 



that is centered at x t and has a width a for each parameter of eq. (jlj). Moreover, we assume 
flat priors and take P(x) = e~ x ^l 2 . A parameter point is hence accepted in the Markov 
Chain, x t+1 = x c , if 

3 (x 2 (z')-X 2 (z c ))/2 



P = 

The x 2 function is computed as 



> 17(0,1). 



Xi{x) - X t 



(6) 



(7) 



where Xi and cr(Xi) are the nominal values of the observables and their l-a errors, 
and Xi(x) are the corresponding values obtained at the parameter point x. Note that, 
including the rr cross section, we have seven measurements (i = 1,...,7) and eight free 
parameters, cf. eq. (jlj). 

We run 5 x 10 5 points on 10 independent chains, with the starting points having very 
different characteristics, such as large or small /i; or M^ 3 smaller, larger, or equal M^. 
We also include two randomly chosen starting points, for which we run 10 6 points. While 
some of the chains start off with huge x 2 ' s ? they all converge fast to x 2 ~ ^(1)- We 
do not only keep the points accepted in the chains but also write out to a separate file 
all points tried which have x 2 < 3.84, corresponding to 95% CL. Figure shows the x 2 
distribution of these points as a function of \i. It has a minimum around /i ~ 600 GeV, 
rises rather steeply for lower values and only slowly for higher values of /i. We can gather 
from this plot that /i is determined only to 500 GeV < /j, < 1200 GeV at la from the 
ILC measurements. The reason is a kind of 'valley' of x 2 minima in the tan /3—fi plane. 
This is illustrated in Fig. Eb which shows in green (black) points which fall within the 2a 
(la) region in the plane of tan (3 versus /i. As can be seen, while tan j3 and \i are not very 
much constrained by the ILC measurements, there is a strong correlation between the 
two parameters. Projections of the la and 2a fitted regions onto the other parameters 
are shown in Fig. [71 The fit results at 95% CL can be summarized as follows. 
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M R3 (GeV) A T (GeV) 

Figure 7: The la (black) and 2a (green) allowed regions projected in different planes: 
(a) M x vs. M 2 , (b) M x vs. fa, (c) M k vs. M L and (d) A T vs. <j> T . 



• The lack of observables individually sensitive to jj and tan (3 leaves a large allowed 
parameter space with 400 < \i < 2200 and 2.8 < tan/5 < 18. Note, however, that 
because the stau mixing is proportional to (A T — /xtan/3) and because in general 
/itan/3 ^> A T these two parameters are strongly correlated. 

• The very precise determination of m^o and m-+ constrains Mi and M 2 to 5 Mi = 
-0.9, +2.1 GeV and 5M 2 = -3.7, +4.9 GeV. 

• Although the phase 0i can take any value, there is a correlation with the value of 
Mi. a large phase is associated with lower values of Mi. This is a direct consequence 
of the precise measurement of the neutralino mass. 

• The measurement of the masses of both staus and their mixing angle well constrain 
Mi 3 and M^ 3 to about 10 GeV. Moreover, solutions with M^ 3 < M^ 3 are strongly 
disfavoured. 
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Figure 8: Predictions for (a) Qh 2 vs. \x and (b) Qh 2 vs. 0i in the la (black) and 2a (green) 
fitted regions. Straight lines show 2a upper limits on Qh 2 = 0.136 



• The trilinear coupling A T is basically undetermined, and the phase is T completely 
unconstrained. This is because the stau mixing is dominated by the term /xtan/3. 
The few points at low tan (3 ~ 4 and /i < 1 TeV in Fig. |Ub have very large 
| Ar | < 3 TeV. Setting an upper limit on \A T \ somewhat tightens the lower bor- 
der of the 2a tan range in Fig. [6b. 

Having determined the allowed parameter space, we next show in Fig.[8]the predictions 
for the relic density of dark matter. The allowed range at 2a is found to be 0.116 < Qh 2 < 
0.19. As expected, Qh 2 tends to be larger than the range that best fits present cosmo logical 
data. The lowest values of Qh 2 are obtained for a large phase <pi and both and vrtf R 
near their maximal allowed value. While there is no strong correlation between /i and 
Mi 2 , see Fig. E](a), it is only for fi w 500-700 GeV that Vth 2 > 0.18. This is because for 
small fi an increased higgsino component (although the LSP is always dominantly bino) 
leads to a smaller annihilation cross section into tau pairs. Furthermore, the mass of 
the LSP is too small for this to be compensated by the very efficient annihilation of the 
higgsino component into W pairs. 

As mentioned above we do not include the light Higgs mass in our fit. This may seem 
surprising since can be measured very precisely at the ILC Moreover, it strongly 
depends on tan /3 and might therefore be used to determine tan f3 and, in turn, constrain 
/i. However, one must not forget that also strongly depends on the parameters of 
the stop sector. Indeed, if the stop parameters are not known, this induces a parametric 
uncertainty of about 15 GeV in nih- As a consequence, the very precise measurement of 
rrih expected at the ILC only poorly constrains tan/5 if the stop sector is not measured 
well. For illustration, Fig. [H](a) shows the correlation between rrih and tan (3 in the 2a fitted 
region for stop parameters fixed to their nominal values. Figure Mjo) shows the same but 
assuming that the stop masses are known to 10% (blue points) or that the stop sector 
is basically unknown (green points). Here note that in our scenario even the discovery 



13 




Figure 9: The 2a region for as a function of tan/3: in (a) for parameters of the stop 
fixed to their nominal values and in (b) for varying stop parameters. In (b), the green 
points show the situation for a completely undetermined stop sector while the blue points 
are for stop masses known to ±10%. 



of the stops is supposed to be very difficult at the LHC because of the overwhelming tt 
background. A measurement of = 116.1 GeV and stop masses known to 10% accuracy 
would constrain tan f3 only to about 3-14, theoretical uncertainties not included. In turn, 
since in our fit the stop parameters are fixed to 1 TeV, we allow Higgs masses as low as 
100 GeV, because one can always find a combination of m^, mj 2 and A t which lifts rrih 
above the LEP limit. 



5 CP-odd observables 

In the previous sections, we have demonstrated how the parameters Mi, M2, fi, etc, of 
our stau-bulk scenario can be determined or at least constrained from an analysis of 
particle masses, cross sections, distributions at the ILC Being CP even, such quantities 
principally allow to determine CPV phases only up to a twofold ambiguity, </> <-» 2ir — <p. 
In our scenario, only a correlation between M\ and <p\ can be obtained, c.f. Fig. 0(b). 
In order to test for non-trivial phases, one needs CP-odd observables. These could be 
T-odd observables based on triple-products in the production and decay of neutralinos or 
charginos, which are kinematically accessible in our scenario, e + e~ — > X1X2 an< ^ e+e ~ ~~ > 
X1X1, respectively! 

5.1 Neutralino production and decay 

For the cross section a of neutralino production, e + e~ —> X1X21 followed by the subsequent 
two-body decay chain X2 ~* f^r T , f x — > Xi T± i we define the T-odd asymmetry [41-43] 

_ a(T > 0) - a(T < 0) 
M ~ a(T>0) + a(T<0y (8) 

2 Note that there cannot be defined any appropriate T-odd asymmetries based on triple products in 
the pair production of staus, since they are scalar particles. 



14 




of the triple product 

T = (Pe- X Pt~) • Pr+, (9) 

of the three- momenta p. In Figure ITDT a). we show the triple-product asymmetry A\ as a 
function of 0i in the la and 2a fitted regions. As can be seen, a measurement of A% ^ 
would be a very clear signal of CP violationjf] Unfortunately, this does not help constrain 
flh 2 due to the two-fold ambiguity in <pi, see Fig. ITOT b). Note also that in our scenario, 
the neutralino production cross sections are suppressed due to the heavy selectron masses 
mg L R = 1 TeV, however can be as large as 10 fb for (P e -,P e +) = (0.8,-0.6) at yfs = 
500 GeV. The branching ratios are of the order of BR(x2 — ► ^t rT ) = 50%. 

In principle, the transverse polarisation of the r from the decay ~^ i s sensitive to 
the phases <f>i, r , and can be large in general [44,45]. However, since the dependence 
of the transverse r polarization on <p T is weak for \A T \ <C tan/3, as in our scenario, we 
only obtain polarisations not larger than 0.7% for non- vanishing phases, with BR(x ] f — > 
f) T T + ) w 50%. 



5.2 Note on chargino production 

CP-odd observables in chargino production and decay are in general only sensitive to the 
phase (p^ of the chargino sector. Our scenario actually assumes M = 0. For completeness, 
we briefly mention how to constrain the phase M in principle. Due to the heavy sneutrino 
mass = 1 TeV, the destructive sneutrino interference term in chargino production 
e + e~ — > XiXi is suppressed, such that the cross section reaches more than 2000 fb in our 
scenario for (P e -,P e +) = (—0.8,0.6) at y/s = 500 GeV. Since all couplings in diagonal 
chargino pair production are real [46], there is no CP-sensitive contribution from 0^ in 
the production process alone. However, if the chargino decays into polarised particles, 

3 The asymmetry A\ flips sign for cf>i £ (tt,2tt), while the CP-even observables, including £lh 2 , are 
symmetric around it. 
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viz. the r, CP-sensitve observables can be denned [47]. Similar to the neutralino decay, 
one such quantity is the the transverse polarisation of the r in the decay Xi ~^ r ^r- In 
our scenario, the transverse r polarisation could be as large as 4% for M = 90°. 

6 Constraining further the model 

In the previous section we have found that the precise determination of a few observables 
at the ILC still leaves a large uncertainty in the prediction of the relic density of the dark 
matter candidate. This holds even if CP-odd observables are included. However, at the 
timescale where the ILC may start its operation, other measurements will be available 
that could constrain further the model. Below we discuss direct and indirect detection 
of DM, as well as EDM measurements. Finally, we comment on further measurements of 
heavy sparticles at the LHC and/or a linear collider with higher energies. 

6.1 Direct detection 

Upper limits on the cross section for elastic scattering of DM particles on nuclei are 
improving every year, with the strongest constraints currently coming from XenonlO [48] 
and CDMS [49], see also [50]. For a 80 GeV WIMP, the limit on the spin- independent 
interaction with protons is a^ 1 < 5 x 10~ 8 pb [49]. The next generation of detectors 
should probe a^ 1 ~ 10 -10 pb within less than 10 years [51]. Furthermore, the near 
maximal sensitivities of the detectors are in the mass range of interest here, i.e. around 
100 GeV. For spin-dependent interactions, limits are much weaker, about 10~ 2 pb [52,53] 
with prospects of reaching ultimately ~ 10~ 6 -10 -7 pb with ton-size detectors, see 
e.g. [54]. 

In our scenario, where the squarks of the first two generations are heavy, the spin- 
independent cross section is completely dominated by the Higgs exchange diagram. We 
therefore expect a^ 1 to be sensitive to /i and <pi through the XiXih coupling. Since the 
LSP couples to the Higgs through its higgsino component, we expect a small cross section 
in our scenario, where the LSP is dominantly bino. This is indeed the case, as can be 
seen in Fig. HH Nevertheless for some of the parameter space the cross section is large 
enough to be detectable at future ton-size detectors. This occurs mostly for near minimal 
fi values, Fig. fTTTa). The dependence on the phase (f>±, Fig. [TlTb). can be related to 
the suppression of the XiX\h coupling. Furthermore, for near maximal phase there is 
an interference between the light and heavy Higgs contributions which leads to a strong 
suppression of the XiP ~> XiP cross section. The models which predict a cross section 
that could be detected in the near future are precisely those where flh 2 is the largest, see 
Fig. [TlT c). Since a positive signal is not expected in this particular scenario, we conclude 
that an improvement on the limit of the neutralino proton elastic scattering would reduce 
the uncertainty in the prediction of the relic density to 0.116 < Qh 2 < 0.17, provided the 
major hadronic uncertainties [55, 56] stemming from the strange content in the nucleon 
can be brought under control. 

The spin-dependent cross section is dominated by Z exchange and hence is also ex- 
pected to be largest when the higgsino component is largest, that is for small /i. No strong 
dependence on the phase <p\ appears. Moreover, as shown in Fig. [TlTd). the prediction for 
the spin-dependent cross section is far below the present limit. The /i-dependence is, how- 
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Figure 12: Indirect detection cross sections av [cm 3 /s] in the la (black) and 2a (green) 
regions, in (a) as a function of /i, and in (b) as a function of 



ever, more pronounced than in the spin-independent case. If ambitious projects like [54] 
can reach ~ 10~ 7 pb, they will indeed probe the small /i region, where our nominal point 
lies. 

Last but not least note that the cross sections for direct DM detection also depend on 
the local DM density and velocity distribution. Therefore even a positive DM signal will 
not directly limit the SUSY parameter space, but rather provide a consistency check. 

6.2 Cross sections for indirect detection 

For completeness we show in Fig. [T2] also the cross section for neutralino annihilation into 
photons. Here we consider only photons that come from the decays of the neutralino 
annihilation products. The cross-section is of the order of few times 10 -26 cm 3 /s and vary 
only by about a factor of two in the 2a parameter range. These cross sections are relevant 
for indirect DM detection experiments such as EGRET [57] or GLAST [58]. While of 
course being interesting and useful by itself, indirect detection will not help in pinning 
down the model parameters. Indeed the rate for 7-rays detection strongly depends on 
the model for the neutralino density near the center of the galaxy, predictions can vary 
by more than an order of magnitude [5] washing out the variations induced by the phase 
dependence. 

6.3 EDMs 

We argued in Sect. [2] that EDM constraints are satisfied in our scenario by choosing heavy 
sfermions of the first and second generations. However, scanning over the parameters of 
the model we find that for not too large values of [i and a large phase 0i, the predictions for 
the Thallium EDM are above the present bound dy\ < 9 x 10~ 25 e cm [59], see Fig. [THlfa.b). 
Taking this bound into account would, however, not much impact the allowed range for 
Qh 2 as can be seen in Fig. [T3l(c). Note that here c?ti is computed using Mg^ ~ 1 TeV and 
that a lack of discovery of sleptons at LHC/ILC will only put a lower bound on the slepton 
mass. Increasing the slepton masses beyond 1 TeV further weakens the EDM constraint. 
For this reason we do not include this contraint in our fit but use it only a posteriori. The 



18 




1000 2000 100 0.1 0.15 0.2 

/i(GeV) Qh 2 

Figure 13: Predictions for d^i in the la (black) and 2a (green) regions, in (a) as a function 
of \i and (b) as a function of <f>±. Frame (c) shows the correlation between g?ti and Qh 2 . 

electron EDM is currently directly extracted from the Thallium measurement, prospects 
for improving the accuracy on the electron EDM by two orders of magnitude are being 
explored [60,61]. Clearly such refined measurements should probe most of the parameter 
space of the scenario considered. 

6.4 Heavy particles at colliders 

A large part of the uncertainty in the model reconstruction comes from the lack of knowl- 
edge of /j, and tan f3. Here the obvious ways out would be to i) determine tan (3 from the 
Higgs sector and/or ii) determine \x from a measurement of the higgsino states. 

As explained above, to exploit the precision measurement of the light Higgs mass, 
one would need to know the parameters of the stop sector. Extracting a stop signal is, 
however, notoriously difficult at the LHC, at least unless t\ is very light [62,63]. This is 
due to the sheer overwhelming ti background. Moreover, even if the stop masses could 
be measured, this would not be sufficient; one would need in addition a measurement of 
the stop mixing angle to extract A t . This only seems feasible at a multi-TeV e + e~ linear 
collider, such as CLIC [64]. 

The heavy higgsino-like chargino xf an d neutralinos X34 could be detected at the 
LHC through Drell-Yan production or through squark decays into them. The former 
suffers from a very small cross section of only few fb, and the latter from tiny branching 
ratios. It is clear that any such measurement will be very challenging and require very 
high luminosity. 

Alternatively, xf or X3 4 could be produced in e + e~ annihilation with higher energy. 
The first kinematically accessible process is X1X34 at about 700 GeV, followed by xtxS 
at about 800 GeV. However, the production cross sections are small, below a femtobarn, 
since one state is almost purely bino or wino and the other higgsino. To be precise, at 
y/s = 800 GeV one would have a{e\ j e~ R — > X1X2) = 0-^4 fb and a{e^e~ R — > X2X3) = 
0.45 fb, increasing to 1.48 fb and 0.74 fb respectively at ^/s = 1 TeV. At ^/s > 1.22 TeV, 
one could finally produce the heavier states in pairs, i.e. xfxfj X3X3, e tc- The X2X2 
cross sections rises fast, giving cr(e^e^ — > X2X2) — fb at y^s > 1.4 TeV, while 
cr ( e fi e Z ~~ * X3X4) — 30 fb at this energy and the other cross sections much smaller. 
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With such measurements, the complete neutralino/chargino system can in principle be 
reconstructed [65-68]. 

A detailed analysis of LHC measurements or linear collider measurements at TeV 
energies is beyond the scope of this paper. Notice, however, that for a percent-level 
collider determination of Qh 2 , one would need to know not only [i and tan/3 but also the 
CP phases (in our case <pi) to a comparable, i.e. percent-level, precision. 

7 Conclusions 

Considering the general MSSM with CP phases, we have investigated a "stau-bulk" sce- 
nario, in which only gauginos and staus are light and accessible to ILC measurements, 
while all other sparticles are heavy. In this scenario, a precise determination of SUSY 
particle properties is quite challenging at colliders because of the dominance of the chan- 
nels involving taus and missing energy. Combining threshold scans, endpoint methods, 
measurements of polarized cross section and tau polarization, the masses of the lightest 
neutralino, chargino and of both staus can be determined precisely at the ILC. Information 
on the stau mixing can also be obtained. 

From these measurements, some of the underlying Lagrangian parameters can be 
extracted with good precision: Mi, M 2 , Mi 3 , M^. Moreover, the product /itan/5 can be 
constrained. However, we are left with a large uncertainty in both /i and tan (3, as well as 
in the phase of Mi, <fti. This causes a rather large uncertainty in the collider prediction 
of the neutralino relic density of 0.116 < Q^h 2 < 0.19 at 95% CL. 

Taking into account the possibility of non-zero phases is of particular importance: 
from pure sparticle spectroscopy, e.g. the measurement of the Xi-fi mass difference which 
is ~ 10 GeV in our scenario, one might conclude that the neutralino relic density is 
too large, Q%h 2 > 0.14. Evidence of a CP-violating signal in EDMs and/or in collider 
measurements would clearly show that phases have to be taken into account. It would, 
however, not directly add information for infering the relic density of the neutralino dark 
matter candidate. 

We have also discussed implications for (in) direct dark matter searches. While the 
cross section for indirect detection shows a too weak dependence, information from large 
scale dark matter detectors could indeed somewhat reduce the allowed parameter range. 
In any case, direct and indirect detection are important cross checks that the x? is indeed 
the DM. 

Finally, for a percent-level collider determination of Qh 2 , which matches the precision 
of cosmological observations, one would need know also /i, tan (3, and the CP phases (in our 
case 0i ) to percent-level precision. To this aim, the above-mentioned ILC measurements 
have to be complemented by precision measurements of the heavy higgsino-like neutralinos 
and charginos at TeV energies. 
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